Prospect for measuring the branching ratio of 
Bs fifi at LHC6 



ELIAS LOPEZ ASAMAR 

Universitat de Barcelona 

Avinguda Diagonal 647, 08028 Barcelona (Spain) 
E-mail: elopez@ecm.ub.es 

ON BEHALF OF THE LHC6 COLLABORATION 

Abstract. The Standard Model predicts a branching ratio for the decay mode Bs fifJ. 
of (3.32±0.32)xl0 while some SUSY models predict enhancements of up to 2 orders of 
magnitude. It is expected that at the end of its life the Tevatron will set an exclusion limit for 
this branching ratio of the order of 10 leaving one order of magnitude to explore. The efficient 
trigger, excellent vertex reconstruction and invariant mass resolution, and muon identification 
of the LHCfe detector makes it well suited to observe a branching ratio in this range in the first 
years of running of the LHC. In this article an overview of the analysis that has been developed 
for the measurement of this branching ratio is presented. The event selection and the statistical 
tools used for the extraction of the branching ratio are discussed. Special emphasis is placed on 
the use of control channels for calibration and normalization in order to make the analysis as 
independent of simulation as possible. Finally, the expected performance in terms of exclusion 
and observation significance are given for a set of values of integrated luminosities. 



1. Introduction 

In the Standard Model (SM), flavour changing neutral currents (FCNC) can be generated only 
through loop diagrams, resulting in low branching ratio (BR) predictions. Some models beyond 
SM can introduce contributions of size comparable to that of SM, yielding significantly different 
predictions for these branching ratios. 

This is the case of the decay Bs — > /x/i, suppressed by helicity, for which the SM predicts 
BR{Bs fifJ.) = (3.32±0.32)xl0 fT. In the Minimal Supersymmetric extension of SM 
(MSSM) this quantity is proportional to tan^/3 [2] , leading to enhancements of up to one order 
of magnitude in cases such as the Non-Universal Higgs Masses framework (NUHM) [3J. 

The current experimental hmit for this BR is BR(i?s fJ-fJ.) < 47x10"^ at 90% C.L. |1] 
(CDF collaboration, with 2 fb^^), roughly a factor 15 above the SM prediction. This channel 
then offer the possibility of observing hints of Physics beyond SM, or in case of confirming the 
SM prediction, rejecting an important region of the parameter space of some of these models. 

2. Experimental conditions 

The LHC 6 detector [5i is specially designed for studying B meson decays produced in LHC from 
pp collisions at a centre of mass energy of ^/s = 14 TeV/c^). The nominal integrated luminosity 
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Figure 1 . Distribution of geometrical likelihood for the signal (red line) and background coming 
from B decays (green, filled), satisfying only the selection cuts (no trigger is applied). 



at the LHC6 interaction point will be 2 fb"^/year, resulting in ~4xl0^^ B meson pairs per year 
inside the acceptance of the detector. The performances most relevant for the measurement 
are the vertexing capabilities to identify the displaced vertices of B decays (impact parameter 
(IP) resolution of 14+35/pT /^m, pT standing for the transverse component of the momentum), 
the invariant mass resolution (~20 MeV/c^, compared to ~35 MeV/c^ at CMS, or ~80 
MeV/c^ at ATLAS), and the muon identification [7J (94% efficient with 1% pollution coming 
from mis-identified n/K). 

3. Analysis strategy 

The first step of the analysis consist of a very loose selection that focus on efficiency for signal 
events, rather than on background rejection. This preliminary selection is based on the quality 
of the vertex, the invariant mass, the impact parameter (IP) of the muons, and the IP/cr/p, 
decay length and momentum of the reconstructed B. It is 65% efficient over reconstructed 
signal events, keeping the bulk of the sensitive signal, and yields ~40 signal events (assuming 
SM predictions) and ~13xl0^ background events per fb~^ under nominal conditions, inside a 
mass window of 60 MeV/c^ around the mass of the Bg. 

The remaining events are caracterized according to three discriminant variables [S]: 

• Invariant mass. 

• Particle-ID likelihood (PIDL) which for both muons combines the likelihood of being muon 
and the likelihood of not being muon. 

• Geometrical likelihood (GL) which combines the information of the distance of closest 
approach of the tracks of the two muons, the minimal IP / a jp of the muons, the isolation of 
both muons, and the lifetime and impact parameter of the -B [S]. Its distribution for signal 
and background is shown in Figure [T} 

Backgrounds are distributed differently from signal in the three-dimensional space spanned 
by the invariant mass and both likelihoods, thus defining regions of different sensitivity to the 
signal. Only backgrounds lying in sensitive regions can affect significantly the measurement. 
Defining the sensitive events as those having GL > 0.5, the background composition of this 
subsample is shown in Table [T] [8 . In this case, the signal-to-background ratio rises from 0.3% 
to 11%. 



Channel 


Yield (2 fb-^) 


{b ^ fi,b ^ fi) 


170 ± 90 


B+ ^ J/'i'Mfi+u^ 


< 20 (90% C.L.) 


B ^ h+h- mis-ID 


8 ± 2 


Signal 


20 ± 2 



Table 1. Sources of background for sensitive events (GL > 0.5). 



4. Calibration and normalization 

Calibration and normalization procedures foreseen for this analysis rely entirely on data, in order 
to minimize dependence on simulations |i8j. Mass sidebands are used to calibrate background 
distributions, while control channels are used for calibration of signal properties: 

• The particle-ID likelihood is calibrated using inclusive samples of J/^^fifj,) (for muon 
hypothesis) and A{pK) (for non-muon hypothesis), which can be selected with high purity 
using only kinematical cuts. 

• The geometrical likelihood and invariant mass are calibrated using B h^h^ , which has 
the same kinematical properties as the signal. The trigger introduces strong biases in the 
distribution of the geometrical likelihood, which can be removed by using events triggered 
on particles not related with the signal (~7% of total B h~^h~). 

Control channels with known branching ratios allow to normalize the signal trough 

iV. = ^f^iVc . (1) 
ec jB,cBnc 

Subindices S and C stand for signal and control channel respectively. The total efficiency e can 
be split as e^f,c/prod ^ ^sei/rec ^ ^trig/seh where erec/prod is the reconstruction efficiency on produced 
events, Ssei/rec is the efficiency of the offiine selection on reconstructed events, and etrig/sei is the 
efficiency of the trigger on offiine-selected events. B^ decays are aproppiate for calibration due 
to the accuracy in the measurements of their branching ratios. A proper choice of the control 
channel leads to cancellation of the effect of some sources of inefficiency: 

• B^ J/^(fifj,)K'^ . The effect of the trigger is similar to that for the signal due to the 
J/^ muons. The ratio of reconstruction efficiencies needs to take into account the different 
number of tracks in the final state between the signal and the control channel (two and 
three tracks respectively), and it is estimated through: 

er-ec(2 tracks) erec(3 tracks) , , 

~ (2) 

ereci'i tracks) erec{^ tracks) 

Using B^ J /'^{iiij)K*^{'kK) as the four-track decay channel the ratio between the left- 
hand side and the right-hand side is 92%. 

• B ^ h^h^ . Has the same kinematic properties to that of signal, leading to cancellation of 
reconstruction and selection effects. 

In both cases the ratio of trigger efficiencies is estimated using events triggered on particles 
not related with signal, which have a relatively small trigger bias on signal properties. 

The ratio of B production fractions in Equation [T] is the main source of uncertainty of 
the measurement, as long as channels coming from i? decays are used for normalization. It 
introduces a systematic error of 13%. 




Figure 2. Expected performance for exclusion (left) and measurement (right) of a given 
branching ratio as a function of the integrated luminosity. The expected limit set by CDF 
at the end of its life (8 fb^^ of data) is also shown. 



5. Extraction of the branching ratio 

The expected distributions of invariant mass and both particle-ID and geometrical likelihoods 
for signal and background will be obtained from the control channels and the mass sidebands 
respectively. These distributions will be combined together and compared with the measured 
ones, using the CL., method [9]. In the absence of signal this analysis would reach the expected 
CDF limit with only an integrated luminosity of 0.1 fb^^, and would reach the SM prediction 
with 2 fb"^ [8] (see Fi gurepj). In case of observing a signal, the 3a {5a) measurement could be 
achieved with 3 fb"^ (10 fb"^) [8 of data. 



6. Conclusions 

The method developed for the measurement of the branching ratio of Bg nfi at LIIC6 makes 
an extensive use of the control channels in order to be independent of simulations. This analysis 
is expected to lead to relevant results even with early data, excluding the SM value at 90% CL 
with only 2 fb~^ in absence of signal, or observing it with only 3 fb~^, thus potentially being 
able to provide one of the first evidences of New Physics from the LHC. 
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